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ABSTRACT 

We  study  the  stability  of  spectral  approximations  to  scalar  hyperbolic  initial-boundary 
value  problems  with  variable  coefficients.  Time  is  discretized  by  explicit  multi-level  or  Runge- 
Kutta  methods  of  order  <  3  (forward  Euler  time  differencing  is  included),  and  we  study 
spatial  discretizations  by  spectral  and  pseudospectral  approximations  associated  with  the 
general  family  of  Jacobi  polynomials.  We  prove  that  these  fully  explicit  spectral  approx¬ 
imations  are  stable  provided  their  time-step,  At,  is  restricted  by  the  CFL-like  condition, 
At  <  Const. TV-2,  where  N  equals  the  spatial  number  of  degrees  of  freedom.  We  give  two  in¬ 
dependent  proofs  of  this  result,  depending  on  two  different  choices  of  appropriate  L2-weighted 
norms.  In  both  approaches,  the  proofs  hinge  on  a  certain  inverse  inequality  interesting  for 
its  own  sake.  Our  result  confirms  the  commonly  held  belief  that  the  above  CFL  stability 
restriction,  which  is  extensively  used  in  practical  implementations,  guarantees  the  stability 
(and  hence  the  convergence)  of  fully-explicit  spectral  approximations  in  the  non-periodic 
case. 
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1.  INTRODUCTION. 


We  are  concerned  here  with  fully-discrete  spectral  and  pseudospectral  approximations  to 
scalar  hyperbolic  equations.  In  this  context,  the  spectral  (and  respectively,  the  pseudospec¬ 
tral)  approximations  consist  of  truncation  (and  respectively,  collocation)  of  W-term  spatial 
expansions,  which  are  expressed  in  terms  of  general  Jacobi  polynomials;  Chebyshev  and  Leg¬ 
endre  expansions  are  the  ones  most  frequently  found  in  practice.  In  this  paper  we  prove  that 
such  IV-terms  approximations  are  stable,  provided  their  time-step,  A t,  fulfills  the  CFL-like 
condition,  At  <  Const. N~2.  To  clarify  the  origin  of  such  CFL-like  conditions  in  our  case, 
we  recall  that  the  Jacobi  polynomials  are  in  fact  the  eigenfunctions  of  second-order  singular 
Sturm-Liouville  problems.  Our  arguments  show  that  the  main  reason  for  the  above  CFL 
limitation  is  the  0(N2)  growth  of  the  corresponding  .W-eigen  values  associated  with  these 
Sturm-Liouville  problems. 

The  paper  is  organized  as  follows.  Section  2  includes  a  brief  summary  on  the  properties 
of  Jacobi  polynomials  (and  their  quadrature  rules)  which  are  used  throughout  the  paper. 
In  Section  3  we  state  our  main  stability  theorems  for  Forward  Euler  time-differencing  and 
( pseudo- )spectral  spatial  differencing,  for  constant  coefficients  equations  with  homogeneous 
boundary  conditions.  Section  4  extends  our  stability  results  to  the  inhomogeneous  case.  In 
Section  5  we  discuss  multi-level  and  Runge-Kutta  time  differencing.  Finally,  in  Section  6  we 
show  how  to  extend  our  results  in  the  presence  of  (positive)  variable  coefficients. 

2.  VERY  SHORT  GUIDE  TO  JACOBI  POLYNOMIALS. 

Jacobi  polynomials,  P^,  are  the  eigenfunctions  of  the  singular  Sturm-Liouville  problem 
(2.1a)  ((1  —  x2)w(x)Pjea (x))  Xilw(x)Pjca'0\x )  =  0,  — 1  <  x  <  1, 

with  corresponding  eigenvalues  A*, 

(2.16)  Afc  =  \k(a,(3)  =  k(k  -I-  a  +  /?  +  1). 

Different  families  of  Jacobi  polynomials  are  associated  with  different  weight  functions  w(x), 

(2.1c)  u;(x)  =  w(x\ a,/?)  =  (1  —  x)a(l  +  x)^,  a,(3>—  1. 

In  the  sequel  we  shall  frequently  use  several  properties  of  the  Jacobi  polynomials.  A  brief 
summary  of  these  properties  is  given  below  (consult  e.g.  [13]).  We  start  with  the  well-known 

PROPERTY  1  (Orthogonality).  We  have 

(2.2)  =  0,  if  k. 
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The  derivatives  of  Jacobi  polynomials  are  also  Jacobi  polynomials.  This  is  evident  from 
the  following  property  which  shows  that  {P£a,/9^  }*>o  are  orthogonal  with  respect  to  the 
weight  (1  —  x2)tu(i)  =  w(x ;  a  +  l,/3  +  1)  and  hence 

(2.3)  P£f  =  Const*,  a,fiP(ka+h0+1\  Const*, a, ^  =  I(Jfe  +  a  +  0  +  2). 

PROPERTY  2  (Orthogonality  of  derivatives).  We  have 

(2.4)  =  0.  3  /  *. 

(2-5)  =  ^ll^’llic). 

Indeed,  (2.4)  and  (2.5)  follow  from  integration  by  parts  of  (2.1)  against  Pja’^(x). 

Let  TTff  denote  the  space  of  algebraic  polynomials  with  degree  <  N.  A  useful  consequence 
of  the  last  two  properties  is  provided  by 

LEMMA  2.1  (Inverse  inequality).  For  all  peics  we  have 

(2-6)  ||p#||(i-*a)w(*)  <  \/*n||p||u>(*)>  pe*jv- 

Here  uj(x)  stands  for  an  arbitrary  w(x]  a,  /3)  weight,  and  Ajv  =  Aj \f(a,(3)  is  the  corresponding 
N-th  eigenvalue. 

Remarks. 

1.  The  inequality  (2.6)  can  be  viewed  as  the  algebraic  analogue  of  the  trigonometric 
inverse  inequality, 

(2.7 a)  ||p'||i,’[— *,*]  <  NlblU’J-w,*],  p  =  any  N  —  trigonometric  polynomial. 

This  should  be  contrasted  with  a  similar  L2-inverse  inequality  for  algebraic  polynomials 
where  there  is  a  loss  of  iV2-factor  for  each  derivative  [3], 

(2.76)  ||p,|k,[-i.il  ^  Const.  7V2||p||£,2(-i,i),  p  =  any  N  —  algebraic  polynomial, 

and  this  estimate,  (2.7b),  is  sharp  in  view  of,  e.g.,  pjv(x)  =  S*LoPfc  *  3\x)-  Thus,  the  use 

of  the  different  weighted  L2-norms  in  the  algebraic  case,  (2.6),  is  essential  in  order  to  retain 
a  loss  of  only  ~  N-factor  for  each  derivative. 
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2.  The  inverse  inequality  (2.6)  can  be  viewed  as  an  Z2-weighted  version  of  Bernstein’s 
inequality 

11(1  -  xa)Sp'(*)|U«[-i.i]  <  ^l|p(*)llfi-[-i.il. 

Standard  interpolation  arguments  between  this  X°°-type  estimate  and  the  L2-type  estimate 
(2.6)  yield  for  q  >  2 


(2.7c)  ||(1  -  xJ)5p'(I)lll.^...„[-i.il  <  •  IW*)IUl(.i.  wI-i,il. 

Similar  weighted  i’-type  estimates  apply  to  higher  derivatives. 


V**N- 


PROOF.  Given  p(x)  in  ir/y,  we  will  use  its  Jacobi  expansion,  p(x)  =  Ylk=oakPka^\x) 
and  p'(x )  =  Ylk=oakPka’^  (x)-  Starting  with  the  left-hand  side  of  (2.6)  and  using  (2.4),  (2.5) 
and  (2.2)  in  this  order,  we  obtain 

‘.0  Y 

k=  0  Jk=0 


1— x3)w(x) 


=  E  Vllld*Wlli(.)  <  MKHS)1 


(lhs)’  =  E“2||P< 


We  note  in  passing  that  Lemma  2.1  can  be  generalized  to  higher  derivatives:  successive 
application  of  (2.6)  with  u(x)  =  w(x;  a,  0)  yields 


(2.8)  ||p(fc)(x)ll(i-x’)M*)  <  II  +  3, 0  +  j)  ■  IIKx)llL(x)»  1**N- 

j=o 

This  leads  us  to  a  ‘natural’  definition  of  non-periodic  Sobolev  spaces  equipped  with  finite 
f/^Ti-n°rm,  where, 


(2.9) 


II*:,. 


fc=o 


With  this  in  mind,  we  now  recover  a  sharp  inverse  inequality  familiar  from  the  trigonometric 
setup 


(2.10)  l|p||»i(.)  ^  Const,  •  JW||pII«<.)>  Const  ~  1  +  -^r,  pei xN. 

In  the  above  discussion  we  can  replace  integrals  by  discrete  summations  in  view  of  the 
well-known 


PROPERTY  3  (Gauss  quadrature  rule).  Let  {qs(x)} n>o  be  a  family  of  ir  ^  -polynomials 
orthogonal  w.r.t.  the  u>(x) -weighted  L2  inner  product.  Let  —1  <  X\  <  x3  . . .  <  x #  <  1  be 
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the  N  zeroes  of  qN(x).  Then  there  exists  positive  discrete  weights,  such  that  for  all 

ptn2N-i  vis  have 


(2.11) 


,1  N 

/  u>(x)p(x)dx  =  “jp(xj), 
J~1  j=i 


Pt*2N- 1- 


Remark.  To  compute  the  Gauss  weights  we  set  p(x)  =  in  (2.11).  Since  p(xj)  =  0,  j  /  k, 

(2.11)  yields 

(2.12)  o)j.  =  — -  [  w(x)-N -  ^  dx,  1  <  k  <  N. 

V  q'N(*k)J- 1  K  Jx-xk  ’  ~  ~ 


PROOF.  We  have  p(x)  =  t(x)qn(x)  +  r(x)  for  some  t(x)  and  r(x)  in  ttn-i-  The  choice 
of  weights  in  (2.12)  guarantees  that  (2.11)  is  valid  for  "Ks-i  -polynomials,  for 


span 


f  \ 

1  x  -  xk  / 


=  7TJV-1. 


Kk<Nf 


This  together  with  our  assumption  that  qtf(x)  is  Z/jj^-orthogonal  to  imply 

r\  ^ 

/  uj(x)p(x)dx  =  /  u;(x)r(x)dx  =  2j  ujr(xj)  -  zJ 
J~1  J~l  j=  l  i=i 


□ 


EXAMPLES. 

1.  Gauss- Jacobi  quadrature  rule.  By  Property  1,  (2.11)  applies  to  with 

w(x)  =  i u(x;a,/3).  Hence  there  exist  {u)j  =  w^(a,  (3)}^  such  that 

rl  " 

(2.13)  /  w(x)p(x)dx  =  2^  Wjp(xj),  for  all  pe-K2N-i- 

J-i  j=l 


Remark.  The  Gauss-Jacobi  quadrature  rule  (2.13)  can  be  used  as  a  highly  accurate  quadra¬ 
ture  rule  for  general  smooth,  not  necessarily  polynomial  functions.  The  error  incurred  in 
such  cases  is  governed  by  [4,  p.  75] 

•i  N 

(2.14)  /  w(x)f(x)  —  Y2w}f(xj)=  Const .  f^N\0),  Const.  >0,  |0|  <  1. 

j=i 
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2.  Gauss-Lobatto-Jacobi  quadrature  rule.  By  Property  2,  (2.11)  applies  to  }at>o 

with  ui(x)  =  (1  —  x2)w(x-,a,  0),  and  therefore,  there  exist  {xbj  =  w^(a  +  1,0  +  l)}jLi  such 
that 

,1  N 

(2.15)  /  (1  —  x2)w(x)r(x)dx  =  Wjr(xj),  for  all  re-K2N-\- 

J~1  j=i 

This  is  in  fact  a  special  case  of  the  Gauss-Lobatto-Jacobi  quadrature  rule  which  is  exact  for 
all  pe-K2N+i-  Indeed,  for  all  such  p’s  we  have  p(x )  =  (1  —  x2)r(x)  +  £{x)  with  r(x)  in  7r2//_i 
and  a  linear  £(x)  =  p(-l)~  +  p(l)Lt*.  By  (2.15) 

/-i  w(x)p(x)dx  =  EjLi  Wjr(xj)  +  /ij  w(x)£(x)  = 

=  E"  i  t +  A  M*)4*)  -  EjLi  =  /  +  //  +  "I- 

Thus,  we  have 

1  =  wjP(xj)>  wf  =  ~^~jwf(a  + 1,0+1) 

j= l  1  xj  1 

and  the  two  expressions,  II  -f  III,  amount  to  a  linear  combination  of  p(— 1)  and  p(l) 

II  +  III  =  w%p(x0)  +  u;^+1p(x^+1),  x0  =  -1  <  xx  <  . . .  <  xN  <  1  =  xN+1. 

Hence,  there  exist  {wj  =  w^(a,  0)}jrJ'o  such  that 


,i  ff+i 

(2.16)  /  w(x)p(x)dx  =  ^2  Wjp(xj),  for  all  pe7T2Ar+i .  □ 

■/"1  j= o 


Finally,  we  shall  need  some  information  on  the  behavior  of  the  collocation  points  which 
appear  on  the  right  of  (2.13)  and  (2.16).  We  have 

PROPERTY  4  ^Distribution  of  zerosl.  If  Xj  =  cos  9j  is  the  j-th  zero  of  Pffl0\x),  then 
[1,  p.  287]  N6j  is  a  corresponding  zero  of  Bessel’s  function,  and  hence 

1  —  x2  =  sin2  9j  ~  ConstJ0iV-2  for  jeJ  ~  {1  <  j  <  j0,  N  —  j0  <  j  <  N}. 

Thus,  the  zeros  of  Pff'0\x)  are  accumulated  within  0(N~2)  neighborhood  of  {  —  1,  +1}.  More 
precise  estimates,  e.g.,  [12,  p.  19]  yield 

1  N 2 


(2.17) 


1  —  Xj  2(1  +  a) 


,  1  <  j  <  N. 
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3.  FORWARD  EULER  WITH  HOMOGENEOUS  BOUNDARY  CONDITIONS. 

We  start  with  the  scalar  constant  coefficient  hyperbolic  equation, 

(3.1)  ut  =  aux,  (a:,  t)e[— 1, 1]  x  [0,oo),  a  >  0, 

which  is  augmented  with  the  homogeneous  condition  at  the  inflow  boundary, 

(3.2)  u(l,t)  =  0,  f>0. 

To  approximate  (3.1),  we  use  the  forward  Euler  time  differencing  on  the  left,  and  either 
spectral  or  pseudospectral  differencing  on  the  right.  Thus,  we  seek  a  temporal  sequence  of 
spatial  -k n- polynomials,  vm  =  t/#(x,  tm  =  mAt),  such  that 

(3.3a)  vN(x,  tm  -f  At)  =  vN(x,  tm)  +  At  ■  v'N(x,  tm)  +  At  •  r(tm)gAr(x). 

Here,  qt/(x)  is  a  ttjv -polynomial  which  characterizes  the  specific  (pseudo)spectral  method  we 
employ,  and  r  =  r(tm)  is  a  free  scalar  multiplier  to  be  determined  by  the  boundary  constraint 

(3.3 b)  v#(x  =  1,  fTO)  =  0. 

We  shall  study  the  spectral-Jacobi  tau  methods,  [8], [2],  where 

(3.4)  vN(x,  tm  +  At)  =  vN(x,  tm)  -f  At  •  av'N(x,  tm)  +  At  ■  T(tm)qN(x),  qN(x)  =  Pjf'0\x), 

and  the  pseudospectral-Jacobi  methods,  [5], [2],  which  are  collocated  at  the  interior  extrema 
pia>P)  •  „ 

OI  r N+i  >  1-e-j 

(3.5)  vN(x,  tm  +  At)  =  vN(x,  tm)  +  At  •  av'N(x,  tm)  +  At  •  T(tm)qN(x),  qN(x)  =  (x). 


Remark.  These  two  families  of  spectral  and  pseudospectral  Jacobi  methods  are  closely  related 
since  Pj^_^  ( x )  is  a  scalar  multiple  of  P^*+1’^+1^(x),  consult  (2.3).  We  will  not  discuss  here  a 
different  alternative  to  (3.5)  where  one  collocates  at  the  interior  extrema  of  Pj^'^(x)  together 
with  the  downstream  outflow  boundary  so  that 

(3.6)  qN(x)  =  (1  +  x)PfrP)  (x). 

Let  —  1  <  xj  <  x2  <  . . .  <  Xff  <  1  be  the  Y-different  zeros  of  the  forcing  polynomial 
qyv(x).  The  spectral  approximation  (3.3a)  restricted  to  these  points  reads 

(3.7a)  vN(x}>  tm+1)  =  vN(xj,  tm)  +  At  ■  av'N(x„  tm),  1  <  j  <  N, 
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and  is  augmented  with  the  homogeneous  boundary  condition 


(3.74) 


v„(l,  ("*)  =  0. 


Equations  (3.7a),  (3.7b)  furnish  a  complete  equivalent  formulation  of  the  spectral  approxi¬ 
mation  (3.3a),  (3.3b).  An  essential  ingredient  in  a  stability  theory  of  such  approximations 
lies  in  the  choice  of  appropriate  £2-weighted  norms 

(3.8)  ||/(x)||2  =< /(i),/(x)  >,  <  f(x),g(x)  >=  '52ujf(xj)g(xj). 

j=i 

We  make 

DEFINITION  3.1  (Stability).  The  approximation  (3.7a),  (3.7b)  is  stable  if  there  exist 
discrete  weights,  {ujj  >  0}^,  and  a  constant  t)0  independent  of  N,  such  that 

(3.9)  ||uw(-,t)||u,  <  Const. eT*’t||ujv(-,0)||w. 

The  approximation  (3.7a),  (3.7b)  is  strongly  stable  if  (3.9)  holds  with  Const.  =  1  and  T]o  <  0, 
i.e.,  if 


(3.10) 


|t>jv(-,<)IU  -  IIuaK‘>o)IL- 


We  recall  that  in  the  Jacobi-type  spectral  approximations  (3.4)  and  (3.5),  the  nodes 
{xjJyLj  are  the  zeros  of  Jacobi  polynomials  associated  with  the  Gauss  and  Gauss-Lobatto 
quadrature  rules.  We  use 


(3.11) 


Ax^  -  min(l  +  xu  1  -  xN ) 


to  measure  the  minimal  gridsize  associated  with  these  Gauss  nodes.  Our  choice  of  discrete 
weights  for  the  stability  of  the  spectral  and  pseudospectral-Jacobi  methods  (3.4), 

(3.5)  will  be  specified  later  on;  these  weights  are  related  (but  not  equal)  to  the  corresponding 
Gauss  weights  indicated  earlier. 

With  this  in  mind  we  have 

THEOREM  3.1  (Stability  of  the  spectral  and  pseudospectral  Jacobi  methods).  Con¬ 


sider  the  spectral  approximations  (3.7a),  (3.7b),  associated  with  the  Jacobi-tau  method  (3.4) 
or  the  pseudospectral-Jacobi  method  (3.5).  There  exists  a  positive  constant  rjo  =  770(0,/?)  >  0 
independent  of  N,  such  that  if  the  following  CFL  condition  holds 


(3.12) 


a  \  <  7?0, 
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then  the  approximation  (3.7a),  (3.7b)  is  strongly  stable  and  the  following  estimate  is  fulfilled 
(3.13) 


Notes. 


1.  The  choice  of  L2- weighted  norms.  Theorem  3.1  deals  with  the  stability  of  both  the 
spectral- tau  methods  associated  with  Jacobi  polynomials  Pj^’^(x),  a,  /3e(  — 1, 0),  and  with 
the  closely  related  pseudospectral  methods  associated  with  P/t-fi  (x)>  a>Pe(~  1)0)-  In  each 
case  we  give  two  different  stability  proofs  which  are  based  on  two  different  choices  of  discrete 
L2-weighted  norms;  these  discrete  weights  {o;J}^.1  are  given  by 


(3.14a) 


UJj  = 


1  +  Xj 
1  —  Xj 


w 


]> 


where  =  Gauss  —  Jacobi  weights  in  (2.13), 


(3.146)  UJj  =  (1  -f  Xj)wj ,  where  {wj}^=1  =  Gauss  —  Lobatto  Jacobi  weights  in  (2.16). 


2.  The  CFL  condition.  The  CFL  condition  (3.12)  places  an  0(N~2)  stability  restriction 
on  the  time-step  At,  and  this  stability  restriction  involves  two  factors.  First,  since  we  expand 
our  solution  in  terms  of  the  eigenfunctions  of  the  Sturm-Liouville  problem  (2.1),  the  CFL 
condition  involves  the  corresponding  ( N  —  l)-th  eigenvalue 


(3.15a)  A„_x  =  <  N2,  a,pe{- 1,1). 

Second,  the  spectral  Jacobi  approximation  (3.7a)  is  collocated  at  Gauss  nodes  which  are 
accummulated  within  0(N~ 2)  neighborhoods  near  the  boundaries,  i.e.,  by  (2.17), 


(3.156) 


1  N2 

Almin  ~  2(1  +  a)' 


Thus,  in  view  of  (3.15a),  (3.15b),  the  CFL  condition  (3.12)  boils  down  to 


(3.16) 


At  •  al  2  <  Const., 


Const.  = 


Vo 


1  +  a 

2  +  a 


>  0. 


3.  The  choice  of  a  stability  norm.  The  stability  statement  asserted  in  Theorem  3.1  is 
formulated  in  terms  of  discrete  semi-norms,  ]|  •  ||w,  which  are  oz-weighted  by  either  (3.14a)  or 
(3.14b).  We  note  that  ||  ■  are  in  fact  well-defined  norms  on  the  space  of  7r/^-polynomials 
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satisfying  the  vanishing  boundary  condition  (3.7b),  i.e.,  corresponding  to  (3.14a)  or  (3.14b) 
we  have 

(3.17a)  \\vif(-,t)\\u,=  [  w(x)\^-v2N(x,t)dx,  vN(l,t)  =  0, 

(3.17  b)  \\vN(-,t)\\u  =  J  ^w(x)(l  +  x)v2N(x,t)dx,  vN(\,t)  =  0. 

Moreover,  in  view  of  (3.15b),  one  may  convert  the  stability  statement  (3.13)  into  the 
usual  Z<2-type  stability  estimate  at  the  expense  of  possible  algebraic  growth  which  reads 

(3.18)  HM-.OIU*)  <  7^—  e_,Ka<i!M-»0)IU(*),  IM',0lll(*)  =  /  w(x)v2N(x,t)dx. 

i  r  Ct  J —  \ 


4.  Exponential  time  decay.  Let  us  integrate  by  parts  the  differential  equation  (3.1) 
against  (1  +  x)u.  Thanks  to  the  homogeneous  boundary  condition  (3.2)  we  find 


(3.19) 

and  therefore, 

(3.20) 


cL  CL 

—  (1  +  x)u2(x,  t)dx  <  —  —  /  (1  +  x)u2(x,  t)dx, 

dt  J-i  2  J- i 


K-.OIIi+x  < e  ia1l«(-,o)||a+*. 


This  estimate  corresponds  to  the  special  case  of  the  stability  statement  (3. 13)  for  the  spectral- 
Legendre  tau  method  (a  =  (3  =  0)  weighted  by  (3.14b).  The  exponential  time  decay  indi¬ 
cated  in  (3.20)  and  more  generally  in  (3.13),  is  due  to  the  special  choice  of  w-weighted 
stability  nomas.  The  weights  in  (3.14a),  (3.14b)  involve  the  essential  factors  1  +  Xj 

or  which  amplify  the  inflow  boundary  values  in  comparison  to  the  outflow  ones.  Since  in 
the  current  homogeneous  case,  vanishing  inflow  data  is  propagating  into  the  domain,  this  re¬ 
sults  in  the  exponential  time  decay  indicated  in  (3.20)  and  likewise  in  the  stability  statement 

(3.13) . 

5.  The  inflow  problem.  A  stability  statement  similar  to  Theorem  3.1  is  valid  in  the  inflow 
case  where  a  <  0.  Assume  that  the  CFL  condition  (3.12)  holds  with  rj0  =  rjo((3,a),  then 

(3.13)  follows  with  discrete  weights  ujj  =  Y^r^Wj  or  u)j  =  (1  —  Xj)wj. 

The  rest  of  this  section  is  devoted  to  the  proof  of  Theorem  3.1  according  to  the  various 
cases  outlined  in  the  four  lemmata  below.  We  start  with 
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LEMMA  3.2  (Stability  of  the  spectral-tau  method).  Consider  the  spectral- Jacobi-tau 
method  (3.4).  Then  TheoremZ.l  holds  with 

(3.21c)  iOj  =  1  +  —Wj,  {wj  =  wf(a,/3)}^=1  =  Gauss  -  Jacobi  weights, 

1  I  i 


'  |(1  +  /?),  a  +  /?<0, 

(3.216)  T]0  =  r]o(a,  (3)  =  •  ot,/3e(- 1,1). 

.  |(1  —  o:),  a  +  /3>  0. 


PROOF.  Squaring  of  (3.7a)  yields 

(3.22)  +2A<  •  a  <  vN{-,tm),v'N(-,tm)  >  +(A t  •  a)2|K(-,  <m)||2  = 

=  IM-,  tm) \\l  +  2A t  ■  al  +  (At  •  a)2II, 

and  we  turn  to  estimate  the  two  expressions,  I  and  II,  on  the  right  of  (3.22). 

First  let  us  note  that  since  the  7r^-polynomial  vn(x,  tm )  vanishes  at  the  inflow  boundary, 
(3.3b),  we  have 

(3.23)  vN(x,tm)  =  (1  —  x)p(x)  for  some  p(x)  =  pjv-i(x)e7r^-i- 


Also,  a  straightforward  computation  shows  that 

(3.24)  (u;(z)]~~)  ~  =  [(/5  ~  +  2)  -  (/3  +  a)x]u;(x)  >  47?0u;(x),  1*1  <  1. 

where  r/0  =  T]0(a,(3)  is  given  in  (3.21b). 

Now,  since  \^vN(x,  tm)v'N(x,  tm)er2^-i,  Gauss  quadrature  rule  (2.11)  implies 


i  =  Ew>7±^vN(x,,ryN(xj,tm)=  f  w(x)\±^vN(x,nv'N(x,r)dx. 

~Z[  1  —  xy  *  x 

We  integrate  by  parts  the  RHS  of  I,  substitute  Vtf(x,tm)  =  (1  —  x)p(x)  from  (3.23),  and  in 
view  of  (3.24)  we  obtain 

(3.25)  /  =  ~\Jx  («>(*)  '  xfp2{x)dx  <  -2Vo\\p\\l[xy 

Next,  let  us  consider  the  second  expression,  II,  on  the  right  of  (3.22).  As  before,  we  substitute 
vN(x,tn)  =  (1  -  x)p(x)  from  (3.23)  and  obtain 

II  =  \\v'N{;tm)\]l  =  -  Xi)v\Xi)  -  p(x^2  ~ 

J=1  1  -  xi 
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<  2Su,j(1  -  x%v\xi))2  +  +  —P\xi)  =  H\  +  IJ2- 

i=l  j=  1  L  X3 

Since  (1  —  x2)(p'(x))2eT2N-2,  the  Gauss  quadrature  rule  (2.13)  followed  by  the  inverse  in¬ 
equality  (2.6)  implies 

III  =  2]C  Vji1  -  x2j)(p'(x j))2  =  2IIp'II(i-*’M*)  <  2^W-i||p|£(x)>  pe*N- 

i=  i 

This  together  with  the  obvious  upper-bound 

II2  =  2H^j73^P2(®i)p2(*i)  <  -7——  ||p|ll(«)i 

:=j  I  ***  j  ^^piin 


gives  us 
(3.26) 


11  -  (2A"-' +  aZ“) 


|_  |2 
IP  !«;(*)• 


Equipped  with  (3.25)  and  (3.26)  we  return  to  (3.22)  to  find 


(3.27)  \\vN(;tm+1)\\l<\\vN(.,i 


*m\||2 


2At  •  a 


2rj0  -  At  a  ^ 


2 


The  CFL  condition  (3.21b)  implies  that  the  squared  brackets  on  the  right  are  nonnegative, 

2 


(3.28) 


2tj0  -  At-  a  (xN-i  +  )  >  770  >  0, 

\  AXndn/J 


and  hence  strong  stability. 

In  fact,  one  more  application  of  the  Gauss  quadrature  rule  yields 


(3.29) 


<„)  =  £&! 


>  EjLi  WiT%-v2N(xjy  tm)  =  ||t>w(-,  fm)|| 
The  inequalities  (3.29),  (3.28)  together  with  (3.27)  imply 
(3.30)  ||M-,r+1)||^  <  (1  -  2rj0At  •  a)|M-,tm)ll?, 

and  the  result  (3.13)  follows. 


2 

U/* 


□ 


Next,  we  take  advantage  of  the  rather  general  setup  we  employed  in  Lemma  3.2.  Specif¬ 
ically,  since  Pfa is  proportional  to  p0*+1-^+1),  consult  (2.3),  we  may  use  Lemma  3.2  with 
T]o(ct,/3 )  replaced  by  rj0(a  +  1,(3  +  1)  to  conclude: 
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LEMMA  3.3  (Stability  of  the  pseudospectral  method).  Consider  the  pseudo  spectrcd- 
Jacobi  method  (3.5).  Then  Theorem  3.1  holds  with 

1  +  Xj 


(3.31a) 


(3.316) 


a >j  =  - -Wj,  |u/j  =  wf  (a  +  1,  (3  +  1)|  _x  =  Gauss  —  Jacobi  wei 


weights, 


a 


Vo  =  T)0(a,  /?)  =  -->  0,  a, /?<—!,  0). 


As  mentioned  before,  alternative  proofs  of  Theorem  3.1  are  possible.  For  example,  fol¬ 
lowing  [6,  Theorem  5.1]  one  may  employ  a  stable  norm  weighted  by  Uj  =  (1  +  Xj)wj  instead 
of  the  u>j  =  weights  used  in  (3.21a),  (3.31a).  We  begin  with 

LEMMA  3.4  (Stability  of  the  spectral-tau  method  revisited).  Consider  the  spectral - 
Jacobi  tau  method  (3.4).  Then  Theorem  3.1  holds  with 


(3.32a) 


u>j  —  (1  +  Xj)wj,  jiUj  =  t u?(a,0)}  ^  =  Gauss  —  Jacobi  weights 


(3.326) 


Vo  =  Vo{a,P)  =  < 


(  a 
2  ’ 


a  +  /3+l  >  0, 


{  |(1  -  /?),  a  +  0  +  1<0, 


a,/3e(- 1,0). 


PROOF.  We  square  (3.7a)  -  this  time  using  the  inner  product  weighted  by  (3.32a), 
l|t'Ar(*.<m+1)lll  =  ||vjyr(-,i)lli  +  2Afa  <  vN(-,tm),  v'N(-,  tn)  >  + 

(3.33)  +(A/-a)2|K(-,<m)||2w  = 

=  ||vat(-,  tm)||2  +  2A t  ■  al  +  (A t  ■  a)2 II, 

and  as  before  we  have  to  estimate  the  two  expressions  on  the  right  of  (3.33). 

The  first  expression,  I,  involves  discrete  summation  of  /(x)  =  (1  +  x)u^(x, 
and  since  /(x)  is  a  7r2)y-polynomial,  its  JV-nodes  Gaussian  sum  is  not  an  exact  integral.  The 
essential  observation  here  is  that  f(2N)  =  Const.  >  0  in  this  case,  and  the  error  estimate 
(2.14)  tells  us  that  the  Gauss  quadrature  rule  is  upper  bounded  by 

N  f 1 

(3.34)  I  =  10,(1  +  Xj)vN(xj,  im)v'N(xj,  tm )  <  /  tu(x)(l  +  x)vN(x,  tn)v'N(x,tm)dx. 

i= i 

Let  us  recall  that  the  homogeneous  inflow  boundary  condition  (3.7b)  implies 

(3.35)  Utf(x,  tm)  =  (1  -  x)p(x),  for  some  p(x)  =  p^_i(x)e7r^_!. 
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Also,  a  straightforward  computation  shows 

(3.36)  (u>(x)(l  +  x))'(l  —  i)  =  [(/3  —  a  +  1)  —  (a  +  P  4-  l)x]u>(x)  >  4 tj0w(x),  |x|  <  1 
where  rj0  =  Tj0(a,j3)  is  given  in  (3.32b). 

We  integrate  by  parts  the  RHS  of  I,  substitute  (3.35),  and  in  view  of  (3.36)  we  obtain 

(3.37)  I  =  1  +  x))'(  1  -  x)2p2(x)dx  <  -277o||p||(1_iMx)- 

Next,  let  us  consider  the  second  expression,  II,  on  the  right  of  (3.33).  As  before,  we  substitute 
Vf/(x,tm)  =  (1  —  x)p(x)  from  (3.35)  and  Gauss  quadrature  yields 

11  =  +  *i)[(l  -  xj)p'(xj)  -  p(Xj)  ]2  = 

(3.38)  =  w(x)(l  -  x2)(l  -  x)(p'(x))2  -  2  /ij  w(x)(l  -  x2)p(x)p'(x)dx 

+  I— i  tn(a:)(l  +  x)p2(x)dx  =  7/j  +  II2  +  II3. 

The  inverse  inequality  (2.6)  with  weight  u>(x)  =  (1  —  x)w(x)  implies 

^  ^N-i||p||(i_i)w(x),  ^n-\  =  +  l,/3), 

and  integration  by  parts  of  II2  gives 

II2  +  II3  =  [(u>(x)(l  —  x2))'  +  u/(x)(l  +  x)\p2(x)dx 

—  2||p2|U(*)  =  2£^:1  WjYT^:P2{xi)  <  a^~IIpII(1-x)u;(x)- 

Consequently  we  have 

(3-39)  II  <  (aw_,  +  ■£-)  M?,..*.,. 

Equipped  with  (3.37)  and  (3.39)  we  return  to  (3.33)  to  find 

(3.40)  IM',(""H)||i  <  ||w(-,i")||i  -  2Af  ■  a  [2,0  -  A  t-a  (\N.t  +  -A— )1  Mj,..,.,., 

and  the  result  follows  along  the  lines  of  Lemma  3.2,  consult  (3.27).  □ 

Lemma  3.4  does  not  cover  the  pseudospectral  Jacobi  methods,  since  by  (2.6)  the  corre¬ 
sponding  Jacobi  parameters  a  +  l,/3  +  1  £  (  —  1,0).  However,  the  proof  of  Lemma  3.4  can 
be  carried  out  in  the  pseudospectral  case  if  we  replace  the  Gauss  quadrature  rule  by  the 
Gauss-Lobatto  one.  We  omit  the  almost  identical  details  (which  are  outlined  for  the  variable 
coefficients  case  in  Theorem  6.2  below)  and  state 
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LEMMA  3.5  (Stability  of  the  pseudospectral  method  revisited).  Consider  the  pseu- 
dospectral  Jacobi  method  (3.5).  Then  Theorem  3.1  holds  with 

(3.41a)  ujj  =  (1  +  Xj)wj,  { Wj  =  Wj(&, P)}jLi  =  Gauss  —  Lobatto  —  Jacobi  weights, 

f  -f,  a  +  /?-fl>0, 

(3.415)  T}0  =  Vo(<*,/3)  =  <  a,/3e(-l,0). 

I  |(l-/3),  a  +  /?+l<0. 

Remark.  The  stability  asserted  in  Lemma  3.5  is  stated  in  terms  of  the  discrete  semi¬ 
norm  ||u/v(-,  t)||f,  =  Wj(  1  +  Xj)vjf(xj,  t )  weighted  by  the  interior  Gauss-Lobatto  weights 

However,  taking  into  account  the  homogeneous  boundary  condition  (3.7b)  and  the 
exactness  of  Gauss-Lobatto  quadrature  for  7t2jv+i  polynomials,  this  amounts  to 

AM- 1  ,i 

IKK-,  Olli  =  Y,  Vyi1  +  xj)vn(xj>  t)=  w(x)(l  +  x)v2n(x ,  t)dx. 

i=o 

4.  FORWARD  EULER  WITH  INHOMOGENEOUS  INITIAL-BOUNDARY 
CONDITIONS. 

We  consider  the  inhomogeneous  scalar  hyperbolic  equation 

(4.1)  ut  =  aux  +  F(x,  t),  (®,<)e[-l,  1]  x  [0,  oo),  a  >  0, 

which  is  augmented  with  inhomogeneous  data  prescribed  at  the  inflow  boundary 

(4.2)  u(M)  =  <?(*),  <>0. 

Using  the  forward  Euler  time  differencing,  the  spectral  approximation  of  (4.1)  reads,  at  the 
Y-interior  zeros  of  qtf(x), 

(4.3a)  vN(xj}  tm+1 )  =  vN(xj,  tm )  +  At  ■  av'N(xj,  tm)  +  AtF(Xj,  tm ),  gN(xj)  =  0, 

and  is  augmented  with  the  boundary  condition 
(4.36)  vN(l,tm)  =  g(tm). 

In  this  section,  we  study  the  stability  of  (4.3),  (4.4)  in  the  two  cases  of 
(4.4a)  Spectral  —  Jacobi  tau  method  :  qs{x )  =  Pjy“’^(x),  a,  (3e(  — 1,0), 
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and  the  closely  related 


(4.45) 


pseudospectral  —  Jacobi  method  :  g#(x)  =  (x),  a)0e(  — 1,0). 


To  deal  with  the  inhomogeneity  of  the  boundary  condition  (4.3b),  we  employ  a  device 
introduced  in  [6,  Section  5].  Namely  we  consider  the  7r //-polynomial 


(4.5) 

If  we  set 

(4.6) 

then  Vjv(x,t)  satisfies  the  inhomogeneous  equation 

(4.7a)  VN(Xj,  fm+1)  =  VN(xjy  tm)  +  At  •  aV'N{xh  tm )  +  A tF(Xj)  fm), 


F(x,t)  =  F(z,t)  +  a^-g(t) 


which  is  now  augmented  by  the  homogeneous  boundary  condition 
(4.76)  VN(l,tm)  =  0. 

Theorem  3.1  together  with  Duhammel’s  principle  provides  us  with  a  priori  estimate  of 
||V//(-)t)|]u,  in  terms  of  the  initial  and  the  inhomogeneous  data,  || V^(-,  0)||w  and  ||F(-,  t)\\u. 
Namely,  if  the  CFL  condition  (3.12)  holds,  then  we  have 

(4.8)  ||V»(->0||w<e-»"|(VJ,(-,0)|(.+  E  A(.e-’»"('-'m)||F(.,(")IU 

0<tm<t 

Since  the  discrete  norm  ||  ■  ||w  is  supported  at  the  interior  zeros  of  qn(x)  where  V//(ij,t)  = 
vN(Xj,  t)  we  conclude 


THEOREM  4.1  (Stability  of  the  spectral  tau  and  pseudospectral  Jacobi  methods).  Con¬ 
sider  the  spectral  approximation  (3.3a),  (3.3b)  associated  with  the  Jacobi-tau  method  (4.4a)  or 
the  pseudo  spectral- Jacobi  method  (4.4b).  There  exists  a  positive  constant  tj0  =  T)0(a,/3 )  >  0 
independent  of  N,  such  that  if  the  following  CFL  condition  holds  (consult  (3.12)), 


(4.9) 


At  •  a  (Xff- 1  +  T - ^  5: 

V  Azmin/ 


Vo 


then  the  approximation  (3.3a),  (3.3b)  satisfies  the  stability  estimate 
(4.10) 

|M-,t)L<e-’"»l||vA,(.,0)[|„+  E  A(  • 


0<tm<t 


l|f(ll")l'"  +  o!i^T)fl,(r)l 
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Theorem  4.1  provides  us  with  a  priori  stability  estimate  in  terms  of  the  initial  data, 
uat(-|0),  the  inhomogeneous  data,  F(-,t),  and  the  boundary  data  g{i).  The  dependence  on 
the  boundary  data  involves  the  factor  of  ,  which  grows  linearily  with  N,  as  shown  by 


LEMMA  4.2.  Let  be  the  discrete  weights  given  by  either  (3.32a)  or  (3.41a). 

Then  there  exists  a  constant  independent  of  N  such  that 


(4.11) 


Mi)l 


<  Const.  N. 


Remark.  The  stability  estimate  (4.10)  together  with  (4.11)  imply 
(4.12) 

IM-.  OIL  <  E  At  e— [||F(-,  +  Const.TV Ai  •  ( C 

0«m<< 

An  inequality  similar  to  (4.12)  is  encountered  in  the  stability  study  of  finite-difference  ap¬ 
proximations  to  mixed  initial-boundary  hyperbolic  systems  [9].  We  note  in  passing  that  the 
stability  estimate  (4.12)  together  with  the  usual  consistency  requirement  guarantee  the  spec¬ 
trally  accurate  convergence  of  the  spectral  approximation;  consult  [7]  for  the  semi-discrete 
case. 


PROOF.  We  consider,  for  example,  the  spectral-tau  method  associated  with  gjv(x)  = 
and  with  discrete  weights  Uj  =  (1  +  Xj)wf(a,  f3).  Using  Gauss  and  Gauss-Lobatto 
quadrature  rules  (2.13)  and  (2.16)  in  this  order,  we  obtain 


and  (4.11)  follows  in  view  of 


^WN+  l(a>  P) 


o(a<0y 

N 


(i) 


p(° 

rN 


,0) 


(1) 


Const.  N2. 


Similar  arguments  (which  are  omitted)  apply  to  the  pseudospectral  approximation  associated 
with  qs(x)  =  Pjy  ^  (x)  and  with  discrete  weights  u>j  =  (1  -f  Xj)w^(a,  (3). 


5.  MULTI-LEVEL  AND  RUNGE-KUTTA  TIME  DIFFERENCING. 

In  the  previous  sections  we  proved  the  stability  of  spectral  approximations  which  are 
combined  with  the  first  order  accurate  forward  Euler  time  differencing.  In  this  section  we 
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extend  our  stability  result  for  certain  second  and  third  order  accurate  multi-level  and  Runge- 
Kutta  time  differencing,  which  were  studied  in  [10],  [11]. 

To  this  end  we  view  our  ir^-approximate  solution  at  time  level  t,v(-,t),  as  an  ( N  + 
l)-dirnensional  column  vector  which  is  uniquely  realized  at  the  Gauss  collocation  nodes 
( v{xi,t ),... 
v(xN,t),v(l,t)). 

The  forward  Euler  time  differencing  with  homogeneous  boundary  conditions  (3.4),  (3.5), 
reads 

(5.1a)  v(tm  +  At)  =  [I  +  At  •  aL]v(tm),  a  >  0, 

where  L  is  an  ( N  -f  1)  x  ( N  +1)  matrix  which  accounts  for  the  spatial  spectral  differencing 
together  with  the  homogeneous  boundary  conditions 

(5.16)  Lv(tm )  =  («'(*!,  O,  •  •  • .  Ax„,  <m),  0). 

Remark.  The  construction  of  a  spectral  differentiation  matrix  L  can  be  accomplished  in  one 
of  two  ways.  One  possibilty  is  carried  out  in  the  physical  space,  by  exact  differentiation 
of  the  unique  -kn- inter polant  which  assumes  the  gridvalues  v(xl,  t), ...,  v(xN,  t),  v(l,  t)  at  the 
corresponding  Gauss  nodes.  This  leads  to  full  ( N  +  1)  x  ( N  +  1)  differentiation  matrices 
L,  which  are  recorded  for  example  in  [2].  Spectral  differencing  in  the  physical  space  is 
then  carried  out  by  a  matrix- vector  multiplication  at  the  expense  of  0{N 2)  operations. 
Alternatively,  spectral  differentiation  can  be  accomplished  in  the  transformed  space.  In  the 
particular  case  of  Chebyshev  (pseudo-)spectral  method,  this  leads  to  a  factorization  of  the 
corresponding  differentiation  matrix  L,  whose  matrix-vector  multiplication  can  be  carried 
out  efficiently  by  FFT  requiring  O(NlogN)  operations,  consult  [8], [2]. 

Theorem  3.1  tells  us  that  if  the  CFL  condition  (3.12)  holds,  i.e.,  if 

(5.2)  A t  ■  a  (\N-i  +  £~— )  <  Vo, 

then  I  +  At  ■  aL  is  bounded  in  the  ca- weighted  induced  operator  norm, 

(5.3)  ||7  +  At  ■  aL\\„  <  e_rj0°'At. 

Let  us  consider  an  (s  +  2)-level  time  differencing  method  of  the  form 

(5.4)  v(tm  +  At)=  i28k[I  +  ckAt  aL}v{tm-k),  ck  >  0,  9k>  0,  £>*  =  L 

fc=0 

In  this  case,  v(tm  -I-  At)  is  given  by  a  convex  combination  of  stable  forward  Euler  differencing 
and  we  conclude 
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COROLLARY  5.1  (Multi-level  time  differencing).  Assume  that  the  following  CFL  con¬ 
dition  holds, 

(5.5)  At  ■  a  (A*.,  +  -!_)  <  nM)'  Ct>0,  U01 

Then  the  spectral  approximation  (5.4)  is  strongly  stable  and  the  following  estimate  holds 

(5.6)  IM-.OIL  <  e-"““||^(-,0)|U)  V .  =  rnin^  >  0. 

k  Cfc 

In  [10],  second  and  third  order  accurate  multi-level  time  differencing  methods  of  the 
positive  type  (5.4)  were  constructed.  They  take  the  particularily  simple  form 

(5.7)  v(tm  +  At)  =  6{I  +  c0At  ■  aL]v(tm)  +  (1  -  6)[I  +  c.At  •  aL]u(tm-*)> 
with  positive  coefficients  given  in  Table  1. 


Second-order  time  differencing 

e 

Co 

C» 

4-level  method  (s  =  2) 

3 

4 

2 

0 

5-level  method  (s  =  3) 

8 

9 

3 

2 

0 

Third-order  time  differencing 

5-level  method  (s  =  3) 

16 

27 

3 

12 

11 

6-level  method  (s  =  4) 

25 

32 

2 

10 

7 

7-level  method  (s  =  5) 

108 

125 

5 

3 

30 

-1Z 

Table  1.  Multi-level  methods. 

Similar  arguments  apply  for  Runge-Kutta  time  differencing  methods.  In  this  case  the 
resulting  positive  type  Runge-Kutta  methods  take  the  form 

(5.8 a)  =  [I  +  At-  aL)v{tm), 

(5.86)  t/*>(r*+1)  =  0kv(tm)  +  (1  -  0*)[/  +  AiaJLjt/<*-1>(r+1),  k  =  2,  •  •  • ,  s, 

(5.8c)  u(tm+1)  =  v(,)(tm+1). 
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We  state 


COROLLARY  5.2  (Runge-Kutta  time  differencing).  Assume  that  the  CFL  condition 
(3.12)  holds.  Then  the  spectral  approximation  (5.8a-c)  with  0  <  0*  <  1  is  strongly  stable  and 
the  stability  estimate  (3.13)  holds. 

In  Table  2  we  quote  the  preferred  second  and  third-order  choices  of  [11]. 


Second  order  time  differencing 

02 

03 

Two-step  modified  Euler  (s  =  2) 

1 

2 

- 

Third  order  time  differencing 

Three-step  method  (s  =  3) 

3 

4 

1 

3 

Table  2.  Runge-Kutta  methods 


Remarks. 

1.  The  above  results  can  be  extended  to  include  nonhomogeneous  data  as  well.  We  omit 
the  details. 

2.  In  [10],  [11],  higher  order  accurate  (>  3)  multi-level  and  Runge-Kutta  time  differ¬ 
encing  schemes  were  constructed.  In  the  present  context  (of  constant  coefficient  spectral 
approximations),  they  amount  to  convex  combinations  of  the  stable  forward  Euler  differ¬ 
encing  I  -f  At  •  ah  and  its  adjoint  I  —  At  ■  aL.  The  above  argument  does  not  cover  these 
cases,  however,  since  in  our  case  the  stability  of  I  ±  At  •  aL  is  measured  by  different  weighted 
norms. 

6.  SCALAR  EQUATIONS  WITH  VARIABLE  COEFFICIENTS. 

We  begin  with 

EPILOGUE.  When  dealing  with  finite-difference  approximations  which  are  locally  supported, 
i.e. ,  finite  difference  schemes  whose  stencil  occupies  a  finite  number  of  neighboring  grid  cells 
each  of  which  is  of  size  Ax,  then  one  encounters  the  hyperbolic  CFL  stability  restriction 

(6.1)  — — |a|  <  Const. 

Ax 
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With  this  in  mind,  it  is  tempting  to  provide  a  heuristic  justification  for  the  stability  of 
spectral  methods,  by  arguing  that  a  CFL  stability  restriction  similar  to  (6.1)  should  hold. 
Namely,  when  Ax  is  replaced  by  the  minimal  gridsize,  Axmjn  =  minj  |xJ+i  —  Xj\  =  0(N~2)) 
then  (6.1)  leads  us  to 

(6.2)  A f  •  |a|iV2  <  Const. 

Although  the  final  conclusion  is  correct  (consult  (3.16)),  it  is  important  to  realize  that 
this  “handwaving”  argument  is  not  well-founded  in  the  case  of  spectral  methods.  Indeed, 
since  the  spectral  stencil  occupies  the  whole  interval  (-1,1),  spectral  methods  do  not  lend 
themselves  to  the  stability  analysis  of  locally  supported  finite-difference  approximations.  Of 
course,  by  the  same  token  this  explains  the  existence  of  unconditionally  stable  fully  implicit 
(and  hence  globally  supported)  finite  difference  approximations. 

As  noted  earlier,  our  stability  proof  (in  Theorem  3.1)  shows  that  the  CFL  condition  (6.2) 
is  related  to  the  following  two  points: 

#1.  The  size  of  the  corresponding  Sturm-Liouville  eigenvalues,  =  0(N2). 

#2.  The  minimal  gridsize,  ^  =  0(N2). 

The  second  point  seems  to  support  the  fact  that  Ax^n  plays  an  essential  role  in  the  CFL 
stability  restriction  for  the  global  spectral  methods,  as  predicted  by  the  local  heuristic  argu¬ 
ment  outlined  above.  To  clarify  this  issue  we  study  in  this  section  the  stability  of  spectral 
approximations  to  scalar  hyperbolic  equations  with  variable  coefficients.  The  principle  raison 
d’etre  which  motivates  our  present  study,  is  to  show  that  our  stability  analysis  in  the  con¬ 
stant  coefficients  case  is  versatile  enough  to  deal  with  certain  variable  coefficients  problems. 
We  begin  with  the  particular  example  introduced  in  [8] 

(6.3)  ut  =  —xux,  (x,  i)e[ —  1, 1]  x  [0,  oo). 

We  shall  show  that  the  CFL  stability  restriction  in  this  case  is  related  to  the  0(N2)- size 
of  the  Sturm-Liouville  eigenvalues  (point  #1  above),  but  otherwise  it  is  independent  of  the 
minimal  gridsize  mentioned  in  point  #2  above. 

Observe  that  (6.3)  requires  no  augmenting  boundary  conditions,  since  both  boundaries, 
x  =  ±1,  are  outflow  ones.  Consequently,  the  various  7T/v-spectral  approximations  of  (6.3) 
with  forward  Euler  time  differencing,  read 

(6.4)  Vtf(x,  tm  +  At)  =  vN(x,  tm )  —  At  ■  xv'N(xy  tm). 

We  have 
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THEOREM  6.1  (Stability).  Assume  that  the  following  CFL  condition  holds 

(6.5)  At\N<  1,  Xif  =  N(N+l). 

Then  the  spectral  approximation  (6.4)  is  strongly  stable  and  the  following  estimate  is  fulfilled 

(6-6)  ||utf(-,t)||1_x2  <  ||uAr(-,0)||1_xi. 

PROOF.  Squaring  (6.4)  yields 

(6.7)  -2At(M-,n,^M^m))1-^+(A03|l^(-,im)||?_xJ 

=  IIM-,mi?-x*+2At-/  +  (At)2-//. 

Integration  by  parts  shows  that  the  first  expression,  I,  is  given  by 

(6.8)  I  =  l-  J  {x{  1  -  x2))'u*  (x,  tm)dx  =  ^||ujv(-,  tm)||?_r>  -  j_x  x*v2n(x,  tm)dx. 

Next,  we  write  xv'N  =  (i%)'-  v^,  and  by  the  inverse  inequality  (2.6)  the  second  expression, 
II,  does  not  exceed 

(6.9)  II  =  ||(xujv(x,  tm)'||i_xj  —  2  J  x2v2N(x,  tm)dx  <  (A#  —  2)  J  x2v2N(x,tm)dx. 
Inserting  (6.8)  and  (6.9)  into  (6.7)  we  end  up  with 

(6.10)  ||M-,*m+1)ll?-*’  <  (l  +  At)-||uw(.,tm)ii2_x>  +  At-[(AN-2)At-2).  j\2v2N(x,t)dx. 
The  CFL  condition  (6.5)  tells  us  that  the  contribution  of  the  second  term  is  negative,  for 

At[(A*  -  2)  A<  -  2]  •  J'  x2v2n(x,  tm)dx  <  — A«||w(-,  <m)||2_l3 

and  the  strong  stability  estimate  (6.6)  now  follows.  □ 

We  now  turn  to  disscuss  scalar  hyperbolic  equations  with  positive  variable  coefficients 

(6.11)  ut  =  a(x)ux,  0  <  a(x)  <  an,  (x,  t)e[— 1, 1]  x  [0,oo), 
which  are  augmented  with  homogeneous  conditions  at  the  inflow  boundary 

(6.12)  u(l,0  =  0. 
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fcr  BY 

We  consider  the  pseudospectral- Jacobi  method  collocated  at  the  TV-interior  zeros  of  (x)- 

Using  forward  Euler  time  differencing,  the  resulting  approximation  reads 

(6.13a)  vN(xj,  tm+l )  =  vN{xj,  tm )  +  At  ■  a(xj)v'N(xj,  tm),  (xj)  =  0, 

together  with  the  boundary  condition 

(6.136)  vu(l,tm)  =  0. 

Arguing  along  the  lines  of  Theorem  3.1,  we  have 

THEOREM  6.2  (Stability  of  the  pseudospectral  Jacobi  method  with  variable  coefficients). 
Consider  the  pseudospectral- Jacobi  approximation  (6.13a),  (6.13b).  There  exists  a  constant 
7/o  =  t/o  (a,  /?), 


(  a 
2  > 


a+/?+  1  >  0, 


(6.14a)  t/o  =  7/0(a,^)  =  | 

such  that  if  the  following  CFL  condition  holds 
(6.146)  At 


{  a  +  (3  +  l  <  0, 


a,/3e(-!,0), 


-f  2  max  Jt&L  I  <  rjo 

1  <7<*  1  -  XjJ 


then  the  approximation  (6.13a),  (6.13b)  is  strongly  stable,  i.e.,  there  exist  discrete  weights 

(6.15a)  =  (1  +  {wj  =  w^(a,Q)\N  =  Gauss  —  Lobatto  weights, 

a{Xj)  (  3  J  J=i 


such  that 


(6.156) 


||vjv('»0II»  ^  ||M',0)IU- 


PROOF.  We  divide  (6.13a)  by  y/a(Xj), 

—r====vtf(xj,  tm+1 )  =  -  ~  -Vff(xj,  tm)  +  At  •  yja(xj)  ■  v'N{x,  tm ), 

\Ja(xj) 

and,  proceeding  as  before,  we  square  both  sides  to  obtain 

IM-><m+1)ll’  =  IM-,*")lli+ 

(6.16)  +2At  <  wW(-,tm),t4(-,*m)  >  +(AOa||a(-K(-,*m)lli 

=  \\vN(;n\\l  +  2At-I  +  (Aty-II. 
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The  first  expression,  I,  involves  discrete  summation  of  the  7T2Ar-polynomial  f(x)  =  (1  4- 
x)vfif(x,tm)v'N(x,tm)  and  since  /(±1)  =  0  (in  view  of  (6.13b)),  the  ( N  +  l)-nodes  Gauss- 
Lobatto  quadrature  rule  yields 
2V+1 

1  =  E  Wji1  +  ij)vN(xj,iTn)v'N(xj,  tm )  = 
i= o 


J  tn(i)(l  +  i)uat(i,  tm)v'N{ x,  tm)dx. 


We  integrate  by  parts  the  RHS  of  I,  substitute  ujv(x,fm)  =  (1  —  x)p(x)  from  (3.35),  and  in 
view  of  (3.36)  we  obtain  as  before  (compare  (3.37)), 

(6-17)  I  < 


The  second  expression,  II,  gives  us 

II  =  EjLi  wXxi)(  1  +  *i)[(l  -  XjWixj)  -  p(^i)]2  < 

(6.18)  <  20^  £*  x  Wj(l  +  xj)2(p'(xj)f  +  2  Ejli  uija(xj)(l  +  xj)p2(xj) 

—  2a00II\  -f  2  •  //2. 

The  inverse  inequality  (2.6)  with  weight  a>(x)  =  (1  —  x)tw(x)  implies 

Hi  =  I|P'II(1-X»)(1-*M*)  ^  ^N-l  l|p||(l— x)ti>(x)>  ^N~i  =  ^jv-i(a  +  1,/?) 

and  the  expression  //2  does  not  exceed 


1  4.  x  .  ^+1  a(T  \ 

<  mix^x,)  — I)  ■  E  w,-(l  -  X,)v2{x})  <  2  •  max^.  IIpII(i-*M*)- 


Ih 

Consequently  we  have 
(6.19) 


//  <  2  (a„\N-x  +  2  •  maxN  j  IIpII(i-xM*)- 


Equipped  with  (6.17)  and  (6.18)  we  return  to  (6.16)  to  find 

(6.20) 


a(xj) 


luAf(-,tm+1)||^  <  ||vAf(-,tm)|!^  -  2A t  \2t]o  -  At  I  a0 0\N-i  +  2  max 

L  V  1  —  Xj  / 

and  (6.18b)  follows  in  view  of  the  CFL  condition  (6.14b). 


I|p||(l-x)ui(x)l 


□ 


Notes. 

1.  The  case  a(x7)  =  a  —  Const.  >  0,  corresponds  to  the  stability  statement  of  Lemma 
3.5.  Similar  stability  statements  with  the  appropriate  weights  which  correspond  to  Lemmata 
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3.2,  3.3,  and  3.4,  namely,  u>j  =  and-  uj  =  (1  +  xi)^(i-)>  bold.  These  statements 

cover  the  stability  of  the  corresponding  spectral  and  pseudospectral  Jacobi  approximations 
with  variable  coefficients. 

2.  We  should  highlight  the  fact  the  stability  assertion  stated  in  Theorem  6.2  depends 
solely  on  the  uniform  bound  of  a(xj )  but  otherwise  is  independent  of  the  smoothness  of  a(x). 

3.  The  proof  of  Theorem  6.2  applies  mutatis  mutandis  to  the  case  of  variable  coefficients 
with  a  =  a(x,f).  If  a(xj,t )  are  C1-functions  in  the  time  variable  then  (6.20)  is  replaced  by 

,  w f 

<  (1  +  Const. At)||t/(-,  tm)||wm,  w”1  =  (1  +xj)— - — 

a{Xj,  t  ) 

and  stability  follows. 

4.  We  conclude  by  noting  that  the  CFL  condition  (6.14b)  depends  on  the  quantity 

maxx<j<7v  rather  than  the  minimal  gridsize,  ^7—,  as  in  the  constant  coefficient  case 

(compare  (3.12)).  This  amplifies  our  introductory  remarks  in  the  beginning  of  this  section, 
which  claim  that  the  0(N~2)  stability  restriction  is  essentially  due  to  the  size  of  the  Sturm- 
Liouville  eigenvalues,  Xjf-i  =  0(N2).  Indeed,  the  other  portion  of  the  CFL  condition 
requiring 

a(x,) 

At  ■  2  max  - - <  t)0 

1  <j<N  1  -  Xj 

guarantees  the  resolution  of  waves  entering  through  the  inflow  boundary  x  =  1.  In  the 
constant  coefficients  case  this  resolution  requires  time-steps  At  of  size  .  However,  when 
the  inflow  boundary  is  almost  characteristic,  i.e.,  when  a(l)  ~  0,  then  the  CFL  condition 
is  essentially  independent  of  Ax,^,  for  (6.21)  boils  down  to  At  •  2a'(l)  <  r/0.  In  purely 
outflow  cases  such  as  (6.3),  the  time-step  is  independent  of  any  resolution  requirement  at 
the  boundaries  and  we  are  left  with  the  CFL  condition  (6.5)  stated  in  Theorem  6.1. 
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the  spatial  number  of  degrees  of  freedom.  -W£~give  two  independent  proofs  of  this 
result,  depending  on  two  different  choices  of  appropriate tl.- -weighted  norms.  In 
both  approaches,  the  proofs  hinge  on  a  certain  inverse  inequality  interesting  for  its 
own  sake.  ^result  confirms  the  commonly  held  belief  that  the  above  CFL  stability 

restriction,  which  is  extensively  used  in  practical  implementations,  guarantees  the 
stability  (and  hence  the  convergence)  of  fully-explicit/spectral  approximations  in 
the  non-periodic  case.  //  n 
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